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ABSTRACT 

The optical [N I] doublet near 5200A is anomalously strong in a variety of emission-line objects. We compute 
a detailed photoionization model and use it to show that pumping by far-ultraviolet (FUV) stellar radiation 
previously posited as a general explanation applies to the Orion Nebula (M42) and its companion M43; but, it is 
unlikely to explain planetary nebulae and supernova remnants. Our models establish that the observed nearly 
constant equivalent width of [N I] with respect to the dust-scattered stellar continuum depends primarily on 
three factors: the FUV to visual-band flux ratio of the stellar population; the optical properties of the dust; and 
the line broadening where the pumping occurs. In contrast, the intensity ratio [N I]/Hj6 depends primarily on 
the FUV to extreme-ultraviolet ratio, which varies strongly with the spectral type of the exciting star. This is 
consistent with the observed difference of a factor of five between M42 and M43, which are excited by an 07 
and B0.5 star respectively. We derive a non-thermal broadening of order 5 km s ' for the [N I] pumping zone 
and show that the broadening mechanism must be different from the large-scale turbulent motions that have 
been suggested to explain the line-widths in this H II region. A mechanism is requked that operates at scales 
of a few astronomical units, which may be driven by thermal instabilities of neutral gas in the range 1000 to 
3000 K. In an appendix, we describe how collisional and radiative processes are treated in the detailed model N 
I atom now included in the Cloudy plasma code. 



1. INTRODUCTION 

The optical emission-line spectrum of a photoionized cloud 
has prominent recombination lines (H i. He i, and He ii) and 
collisionally excited lines (forbidden lines such as [O iii], [O ii], 
[N ii] and [S ii]). The forbidden lines are produced by ions that 
exist within the H^ region, where the gas kinetic temperature is 
high enough (~ 10"* K) for the lines to be collisionally excited 
(Osterbrock & Ferland 2006, hereafter AGN3). Ions with 
potentials smaller than H" exist mainly in the photodissociation 
region (PDR), a cold {T < 10^ K) region beyond the H+-H° 
ionization front which are shielded from ionizing radiation. 
The PDR does not produce strong optical emission due to its 
low temperature. 

The [Ni] doublet at AA5\99 A is an interesting exception 
to this rule. Atomic nitrogen has an ionization potential only 
slightly larger than that of hydrogen, 14.5 eV for N", as op- 
posed to 13.6 eV for H" (Moore 1993). These, together with 
the relatively slow charge exchange reactions between H and 
N (Kingdon Sl Ferland 1996), mean that little N" is present in 
warm gas, so [N i] has a small collisional contribution and the 
lines are generally weak. This expectation appears to be con- 
firmed in high-resolution observations of nearby H ii regions 
such as Orion (Baldwin et al. 2000), hereafter B2000, where 
the doublet has an observed intensity of only 3 x 10"^ that of 
HyS. But we show in this paper that the ratio becomes higher 
within the central parts of the Orion Nebula, and much larger 
in the nearby M43 nebula. 

This study is motivated by the exceptionally strong intensity 
of the [N i] doublet in several unusual classes of nebulae. Fil- 
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aments in cool-core clusters of galaxies and filaments in the 
Crab Nebula can have the [N i] doublet nearly as strong as Hy6 
(Ferland et al. 2009; Davidson & Fesen 1985). The great [Ni] 
strength is the single most exceptional spectroscopic feature 
in the optical region for these nebulae, and could indicate that 
atomic gas has been heated to temperatures warm enough to 
collisionally excite the line. This could be done by a large 
flux of very hard photons or energetic particles, but is an area 
of active investigation. Large-scale velocity variations within 
these objects could also enhance the absorption of continuum 
photons, making continuum fluorescence more important. The 
fact that several very difi'erent physical processes may be active 
makes it diflicult to understand what the strong [N i] doublet 
tells us about these unusual environments. It is, therefore, 
important to quantitatively explain these lines in the arguably 
simplest case, an H II region. 

Continuum fluorescent excitation has been proposed to be 
an important contributor to the intensity of the [N i] doublet 
(Bautista 1999). This process is unusual because the ground 
term of N*^ is not connected to the upper levels of the observed 
[N i] doublet by any LS-allowed transitions. It is the break- 
down of LS coupling in N i which makes the process fast. The 
FUV lines which pump the upper levels of the [N i] doublet 
lie in the wavelength range AA95\A - 1 161 A. The resulting 
intensity of optical [N i] lines will depend on the atomic tran- 
sition probabilities (a difficult atomic physics problem due to 
the breakdown of LS coupling), the Spectral Energy Distri- 
bution (SED) of the incident stellar radiation field around the 
AA95\A — 1 161 A driving lines, and gas motions in the region 
where continuum fluorescence occurs since the driving lines 
become self shielded. Appendix A.2 describes the fluorescence 
mechanism in detail. 

The purpose of this paper is to use the Orion star-forming 
region to check whether photoionization simulations can self 
consistently account for the observed [N i] intensity. Orion 
is a relatively quiescent environment that can serve as a test 
bed for conventional nebular theory. Our simulations largely 
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confirm the prediction by Bautista (1999) that the [N i] lines 
are predominantly formed by continuum fluorescent excitation. 
We show that their intensity relative to HjS is mainly set by the 
non-thermal component of line broadening in shallow regions 
of the PDR. The line broadening needed to account for the 
observed line intensities is consistent with that seen in Orion. 

The 5198, 5200A pair of lines are denoted as ^5199+ in 
this paper. These lines often appear as a single feature at low 
resolution or when the intrinsic line widths are large. Ap- 
pendix A. 3 describes how the two lines within /15 199-1- can 
be used to measure density if the [N i] lines are collisionally 
excited. 

2. OBSERVATIONS 

In the Bautista (1999) study of [N i] emission there was only 
a limited attempt to compare the results with observed line 
intensity ratios. This was done in a qualitative way for numer- 
ous planetary nebulae, supernova remnants, and Herbig-Haro 
objects, and it should be noted that the axes in his Figure 4 are 
all 100 times too large. The Orion Nebula (M 42, NGC 1976) 
presents an excellent opportunity for testing theories of the 
formation of [N i] emission as the lines are known to be present 
under various conditions. 

Fortunately, there is a recently published spectrophotometric 
study (O'Dell & Haiiis 2010, henceforth OHIO) covering aU 
of the brightest part of the Orion Nebula (the Huygens Re- 
gion), the fainter outer region (the Extended Orion Nebula), 
and the nearby Hii region M 43 (NGC 1982). The inclusion 
of M 43 is particularly important since that object lies along 
the borderline between an object being a photoionized Hii 
region and its being a simple reflection nebula. This status 
is caused by the dominant star NU Ori (spectral type BO. 5, 
O'Dell et al. 201 1) being much cooler than the dominant ion- 
izing star of M 42 (6^ Ori C, spectral type 07 V, O'Dell et al. 
201 1). OHIO obtained moderate spectral resolution long slit 
samples at various distances from fl' Ori C and NU Ori. Red- 
dening corrections were determined for each spectrum. In 
additon to emission-line ratios relative to Hfi, absolute surface 
brightnesses in H/3, were determined. An important measure- 
ment made in OHIO was that of the underlying continuum, the 
strength of this continuum being expressed as the equivalent 
width (EW(H^) = /(H/3)//(Cont), where I(U/3) is the surface 
brightness in the H/3 emission-line and /(Cont) is the surface 
brightness of the observed continuum per Angstrom). The 
units for EW(HyS) are Angstroms (A). The expected EW(Hy8) 
for the Huygens Region due to atomic processes is about 
1700 A (O'Dell 2001). It has long been known (Baldwin 
et al. 1991) that the observed equivalent width (EW(I^, Obs)) 
is much smaller than this. This indicates a strong scattered 
light component arises from Trapezium starlight backscattered 
by dust lying in the dense photon dominated region (PDR) that 
lies just beyond the ionized layer that separates 0' Ori C and the 
background Orion Molecular Cloud. OHIO demonstrate that 
EW(Hj0, Obs) decreases with increasing distance from 0' Ori C 
and that EW(HyS, Obs) values for M 43 are comparable to the 
more distant samples within M 42. OHIO determined that the 
M 42 spectra beyond about 1 0' are increasingly affected by 
scattered light originating from the Huygens region. We have 
included only those samples from their "inner" region group 
with distances of less than 8' and all of their M 43 samples in 
this analysis. 

The reddening corrected results from OHIO are shown in 
Figure 1 . Panel A presents the reddening-corrected emission- 



line ratio /([Ni])//(Hy6) as a function of distance from the 
dominant star (0' Ori C for the M 42 results and NU Ori for 
the M 43 results), where /([Ni]) is the total emission from 
the forbidden N^ lines near 5200 A. Panel B presents the 
ratio /([N i])//([0 1]) (where /([O i]) is the sum of the neutral 
oxygen lines at 6300 A and 6363 A) as a function of dis- 
tance. Panel C presents EW(H^, Obs) as a function of distance. 
Panel D presents in logarithmic scale the /([Ni])//(HyS) ratio 
as a function of EW(Hy8, Coit), where EW(Hj6,Coit) is the 
EW(HyS, Obs) value corrected for the expected atomic contin- 
uum of 1700 A. 

In Figure 1 , panel A we note that although there is a wide 
scatter, there is a general increase in /([Ni])//(Hy6) with in- 
creasing distance from 0' Ori C. The M 43 line ratios are much 
larger and show an even more rapid increase with distance from 
NU Ori. In panel B we note a small general increase in the 
^([Ni])/^([Oi]) ratio with increasing distance from 0' Ori C, 
while the three samples for M 43 show ratios much larger than 
those for M 42. There are fewer samples for /([O i]) in M 43 
because of its much lower surface brightness. The most distant 
ratios in M 42 show a large scatter because of the difficulty in 
separating faint nebular emission from the strong foreground 
night-sky [Oi] emission. Panel C shows that EW(Hj0, Obs) 
decreases markedly in M 42 (the scattered light continuum 
becomes stronger with increasing distance from 0' Ori C). The 
scattered light continuum is always stronger in the M 43 sam- 
ples, but there is no obvious correlation with distance from 
NU Ori. We should note here that the blister model for M 42 
is well established, so that we can expect a monotonic change 
in conditions when looking at lines of sight of greater distance. 
However, the physical model for M 43 is not established. Is 
M 43 the simple Stromgren sphere with overlying foreground 
material in the east as suggested by its circular appearance or 
is it too a blister model object? 

In Figure 1, panel D we see that the the M 42 and M 43 
samples form a well-defined sequence when considering the 
/([N i])//(H/3) versus EW(Hy6, Coit). A linear relation would 
indicate that the ratio /([N i])//(Cont, Corr) (which we will 
call EW([N i], Corr)) is constant. /(Cont, Corr) is the observed 
continuum corrected for the atomic component. Considering 
the two nebulae separately, we calculate EW([Ni],Corr) to 
be 1.98 ± 0.65 for M 42 and 2.15 + 0.99 for M 43. The 
presence of an approximate linear correlation suggests that 
[N i] emission is driven by non-ionizing continuum radiation. 
The value of EW( [N i] , Corr) can become a quantitative test 
for any suggested driving mechanism for the [N i] emission 
and is pursued in the remainder of this paper. 

3. PREDICTED EMISSION FROM A RAY THROUGH 
INNER REGIONS OF THE ORION NEBULA 

This highest S/N observations are for bright inner regions of 
the Orion Nebula. To quantify the various physical contribu- 
tors to the formation of [N i] lines we recomputed the (Baldwin 
et al. 1991) model of a ray through the the Hii region. Ap- 
pendix A describes recent improvements in the treatment of N i 
emission in the spectral simulation code Cloudy which we use 
to compute the spectrum. The model is a layer in hydrostatic 
equilibrium: the outward stellar radiation pressure, largely due 
to grains, is balanced by gas, turbulent, and magnetic pressures 
within the nebula. The parameters are those given in BFM 
with the following exceptions: 

• We include the five high-mass stars of the Trapezium 
(see Table 1), using atmospheres from Lanz & Hubeny 
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Figure 1. These four panels present the spectrophotometric results for the sample regions of M 42 and M 43 as described in the text. Filled circles represent M 42 
samples and filled squares represent M 43 samples. The distances are from the center of the Trapezium for the filled circles representing M42 and from NU Ori for 
the filled squares representing M43. 



(2003) and Lanz & Hubeny (2007). This produces more 
1000 A photons relative to the Lyman continuum than 
would be obtained from 0' Ori C alone. 

• We continue the calculation into the PDR and H2 re- 
gion, including the full H2 model described by Shaw 
et al. (2005) and the chemistry network described by 
Abel et al. (2005). The calculation stops at a thickness 
corresponding to Ay = 10^. 

• We work in terms of stellar luminosities and the physical 
size of the blister As a result the model is not plane par- 
allel, it has a ratio of outer to inner radius of about two. 
We simulate observing this structure by using the option 
to integrate intensities along a pencil beam through the 
geometry. 

• The gas is assumed to be in hydrostatic equilibrium, as in 
BFM. We include magnetic, but not turbulent, pressure 
in the gas equation of state. 

• A "tangled" magnetic field is assumed, as described 
in Appendix C of Henney et al. (2005b), with an ef- 
fective magnetic adiabatic index of ymag - 10. The 
magnetic field in the ionized gas is chosen so as to give 
a ratio of gas pressure to magnetic pressure (plasma 



y6) of 10, which is a typical value found for Hii re- 
gions (Heiles et al. 1981; Harvey-Smith et al. 2011; 
Rodriguez et al. 2011). Together with the assumption 
of Tmag = 1 -0, this implies a constant Alfven speed of 
fA - 3.5 km s"', which is roughly consistent with both 
numerical simulations (Arthur et al. 2011) and obser- 
vational limits (Crutcher et al. 2010). The magnetic 
pressure and gas pressure are therefore roughly equal in 
the PDR {T ~ 1000 K,/3~ 1), whereas magnetic pres- 
sure dominates in the colder molecular gas {T ~ 100 K, 
/3« 1). 

Table 1 lists the stars we include. Figure 2 compares two 
SEDs. The lower curve is 0' Ori C by itself while the higher 
curve includes all stars. The largest differences are in the in- 
tensity of the FUV relative to the Lyman continuum. In the 
case where [N i] is photoexcited and Hf) produced by recom- 
bination, the line intensity ratio is proportional to the ratio of 
the FUV relative to the Lyman continuum. The [N i] pumping 
rate will depend on the intensity of the stellar radiation field at 
the wavelengths of the FUV N i lines. Photospheric absorption 
lines are present across the FUV, making an accurate stellar 
model essential. 

The true atmosphere of 0' Ori C remains highly uncertain. 
The object is a close binary with an extended atmosphere and a 
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Table 1 

Massive stars in M42 and M43 



Star 


M/Mq 


SP type 


logL/Lo 


T/K 


log g 


Refs 






M42 inner 








e' Ori A 


14 


B0.5 V 


4.45 


30,000 


4.0 


1 


e' Ori B 


7 


B3 V 


3.25 


18,000 


4.1 


2 


e' Ori C 


32 


07 V 


5.31 


39,000 


4.1 


1, 3 


e' Ori C2 


12 


Bl IV 


4.20 


25,000 


3.9 


3,4 


e' Ori D 


18 


B0.5 V 


4.47 


32,000 


4.2 


1 






M42 outer 








Ori A 


30 


09 V 


4.93 


35,000 


4.0 


1 


Ori B 


7 


B0.5 V 


4.11 


29,000 


4.1 


1 


Ori C 


6 


B4 V 


3.00 


17,000 


4.1 


5,6 


LP Ori 


10 


B1.5 V 


3.75 


23,000 


4.1 


5,6 


PI 744 


5 


B5 V 


2.70 


16,000 


4.1 


5,6 








M43 








NU Ori 


18 


B0.5 V 


4.42 


31,000 


4.2 


7 



Note. — Stellar parameters of all stars more massive than 5 Mq within the 
confines of M42 and M43, divided into three groups. The "M42 inner" group 
are the Trapezium stars, which excite the bright Huygens region of the Orion 
Nebula. The "M42 outer" group are situated 2-10' south of the Trapezium 
and contribute to the excitation of the Extended Orion Nebula. 



References. — 1. Simon-Di'az et al. (2006); 2. Weigelt et al. (1999); 
3. Schertl et al. (2003); 4. Lehmann et al. (2010); 5. Malkov (1992); 6. Fitz- 
patrick & Massa (2005); 7. Simon-Di'az et al. (201 1). 



detected and periodically variable magnetic field. In addition to 
periodic variations with a period of 15.4 days there are known 
non-periodic radial velocity and spectral variations. These 
characteristics are summarized in Stahl et al. (2008). The es- 
tablished complexity of the atmosphere means that predictions 
of the SED of simple atmosphere models have a corresponding 
uncertainty of undefined magnitude. 
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Figure 2. The SED of 6' Ori C is the lower curve while the heavier higher 
curve gives the SED of the Trapezium stars, using the stellar parameters 
summeiized in Table 1 and the predictions of Lanz & Hubeny (2003) and 
Lanz & Hubeny (2007). The figure is centered on 0.1 pm, which is lOOOA. 



3. 1 . Properties of the cloud 

The upper panel of Figure 3 shows the temperature structure 
of the cloud along our ray. The region has a temperature of 
around IC* K while the gas kinetic temperature falls to around 
300 K in the H" region or PDR. The deeper H2 region is also 
colder. The region is thicker than was found in Baldwin 
et al. (1991) due to magnetic support. 

The lower panel of Figure 3 shows the volume emissivity of 
the /15199^ lines along this ray. For reference the lower panel 
also shows the emissivity of some well-observed H2, CO, and 
[O i] lines. We see that both [O i] and [N i] lines form near the 
H^- H° ionization front, while the H2 and CO lines form near 
the H" - H2 dissociation front. 

The H2 line is mainly formed by continuum fluorescent 
excitation for a PDR near an H 11 region (Tielens & HoUenbach 
1985). This formation processes is very similar to that forming 
the [N i] lines. The H2 and [N i] lines form at either edge of the 
PDR due to a combination of abundance and FUV line optical 
depth eff'ects. 
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Figure 3. The temperature, extinction, and emissivity of several important 
emission lines, are shown as a function of depth into the layer. The upper 
panel shows the log of the gas kinetic temperature and the visual extinction 
Ay, with the values of both indicated on the left axis. The lower panel shows 
the normalized emissivity for several important lines. This is the volume 
emissivity (erg cm"' s" ' ) divided by the peak emissivity for each line to place 
them on the same scale. 

Figure 4 shows the volume emissivity of the [N i] line as 
a function of gas kinetic temperature. This is a convenient 
way to visualize the rapid changes in emissivity that occur 
near the H^-H" ionization front, where both emissivity and 
temperature change rapidly. This does not indicate the total 
contribution of various processes to the observed line since 
the surface brightness is the integral of the emissivity over the 
emitting volume. The size of each volume element changes 
dramatically as the conditions change. 
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The peak emissivity occurs at a gas kinetic temperature of 
~ 4000 K, with contributions from continuum pumping and 
collisions. The emissivity increases as the temperature de- 
creases due to the increasing N" abundance in cooler regions. 
Continuum pumping increases with increasing N" abundance. 
The emissivity decreases at lower temperatures due to increas- 
ing optical depths in the FUV lines. The rise in emissivity 
at temperatures ~ 300K is due to formation by molecular 
dissociation, using the estimates outlined in Appendix A.4. 
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Figure 4. The volume emissivity of [Nl] A5199^ vs gas kinetic temperature. 
This is a different view of the data in Figure 3. The total emissivity, and 
the contributions from collisions, continuum pumping, and the chemical 
dissociation processes described in Appendix A, are shown. 



3.2. The effects of non- thermal broadening 

Given these assumptions the only free parameter is the tur- 
bulent contribution to the line width. Figure 5 shows the 
predicted [Ni] 5199^ / Hy6 intensity ratio as function of the 
turbulence. We have determined the broadening of the [N i] 
lines from spectra made with the HIRES spectrograph as part 
of a program studying mass loss from proplyds in the Huygens 
Region (Henney & O'Dell 1999). In the proplyd-free portions 
of these long-slit spectra, the average observed FWHM was 
13.65 + 1.91 km s"'. The instrumental FWHM of the compari- 
son lines was 8.28 ± 0.40 km s"'. If the [Ni] emission arises 
from the region with T — 1000 to 3000 K, then the thermal 
component of the broadening would be 2 to 3 km s After 
quadratic subtraction of the instrumental and thermal widths 
from the observed FWHM, there is a residual non-thermal 
broadening component of 10.6 +1.9 km s The four pro- 
plyds in the sample (150-353, 170-337, 177-341, 182-413) 
have an average distance from 0' Ori C of 0.56 + 0.29'. After 
examination of Panel A of Figure 1 , we see that the expected 
ratio /([Ni])//(Hy8) would be 0.0034 + 0.0007. These values 
are indicated in Figure 5 for comparison with the predictions 
of our model. The linewidth is given as an upper limit, since 
the observed emission profile will in general include contribu- 
tions from macroscopic and microscopic broadening processes 



(see discussion in Appendix C) whereas only the latter will 
contribute to the pumping efficiency. 

In Appendix A we show that in the fluorescent scenario the 
strength of the [N i] emission lines with respect to H/3 is close 
to linearly proportional to the degree of line broadening in the 
region in which the lines are pumped. In order to reproduce 
the observed brightness, our model requires a FWHM for the 
broadening (assuming a Gaussian line profile) of ^ 10 km s"'. 
If this broadening were to be thermal, then a temperature of 
> 10, 000 K in the pumping region would be required, which is 
much larger than the ^ 2000 K predicted by our Cloudy model. 
Instead, it is likely that the majority of the broadening is non- 
thermal in nature. Significant non-thermal line widths have 
been reported in the spectra of Orion Nebula emission-lines 
(O'Dell 2001; O'Dell et al. 2003; Garci'a-Di'az et al. 2008). 
The nature of the processes producing this broadening is not 
known, but it must be important as its magnitude indicates 
that as much energy is contained there as is contained in the 
components explained by basic photo-ionization physics. 
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Figure 5. Surface brightness of the [N I] line relative to H/S as a function of the 
FWHM of the non-thermal line broadening component in the [N I] formation 
region. The observed ratio is indicated along with en'or bars that represent the 
scatter in the ratio. The observed FWHM is indicated and is really an upper 
limit to the microturbulent broadening, as discussed in Appendix B. 



4. DISCUSSION 

4.1. Spatial variation of intensity ratios 

In § 3 we established a model of an essentially substellar 
point in the Orion Nebula that adequately explains the observed 
^([N i])//(Hy6) line ratios. However, we see in Figure 1, Panel A 
that this ratio varies across the Huygens Region and its near 
vicinity, rising monotonically with increasing distance from 
0' Ori C. A similar increase is seen for the line ratios in M43 
with increasing distance from NU Ori. On the other hand. 
Figure 6 shows that the corrected equivalent width of the [N i] 
lines is essentially constant at 2 ± 1 A and shows no detectable 
variation either within M42, nor between M42 and M43. Since 
the PDR is optically thick to the irradiating stellar continuum. 
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its visual scattered light is a measure of the FUV continuum 
that pumps the upper states of 5199+. The quantitative relation 
is determined by the scattering properties of the grains. These 
are elaborated in Appendix B. 




200 400 600 

EW (Hp.Corr) (Angstroms) 

Figure 6. Equivalent width of tlie [N l] 5 199+ versus equivalent width of Ufi. 



Table 2 

Spectral energy distributions over selected wavelength intervals 



Group 



SED ratio 

FUV/visual EUV/visual FUV/EUV EW(H/3, Int), A 



M42 inner 


20.34 


7.36 


2.76 


380 


M42 outer 


18.18 


4.39 


4.14 


213 


M43 


19.02 


1.40 


13.62 


81 


Ring Nebula 


76.33 


192.24 


0.40 




Crab Nebula 


1.10 


1.12 


0.99 





Note. — Columns 2—4 show the ratio of (AL^) between different wave- 
length bands: EUV = 507-912 A; FUV = 950-1200 A; visual = 4800^900 A. 
Column 5 shows the "intrinsic" equivalent width Lh^ / of the stellar pop- 
ulation, where the continuum luminosity L,i is evaluated adjacent to the HyS 
line and line luminosity L^g is calculated as described in the text. Results are 
shown for the three stellar groupings of Table 1, using atmosphere models from 
Lanz & Hubeny (2003, 2007). The Ring and Crab Nebulae are considered in 
§4.5. 

5 km s"' if the illumination and viewing angles are both close 
to face-on (see Fig. 1 lb). The dust scattering effective albedo 
is therefore constrained to be 0.5 ± 0.3, which is consistent 
with the expectations for back-scattering from the background 
molecular cloud (see Appendix B.2.3 and Fig. 12) so long as 
the single-scattering albedo is relatively high. 



Appendix B shows that the line intensity ratios and equiva- 
lent widths can be expressed in terms of the ratios of "scatter- 
ing" efficiencies (albedos) and the ratios of the stellar contin- 
uum luminosities in different wavelength bands. The pumping 
contribution to 5199-1- depends on the intensity of the SEE) 
around 1000 A while the intensity of 11/3, which forms by re- 
combination, is proportional to the continuum intensity at 
hydrogen-ionizing energies. These are denoted by FUV and 
EUV in Table 2. In the simplest case we expect the 5199-i-/Hy6 
intensity ratio to scale with FUV/EUV, while the 5199-1- equiv- 
alent width should scale with FUV/visual. 

4.2. The importance of the FUV / EUV ratio 

Table 2 lists the ratios of the average value of the spectral 
energy distribution (SED) AL^i, calculated for the visual, FUV, 
and ionizing EUV bands, and for three different OB stellar pop- 
ulations characteristic of the inner Orion nebula (Trapezium 
region), the outer Orion nebula, and M43 (we will consider 
the Crab and Ring Nebulae, objects very different from Orion 
and with strong [Ni] emission, in § 4.5 below). It can be seen 
from the table that the FUV/visual luminosity ratio is approxi- 
mately constant between the three stellar populations in Orion 
(variation < 10%), which arises because the spectral shape in 
all cases approximately follows the Rayleigh- Jeans behavior 
of ~ A . On the other hand, the relative strength of the 
EUV band with respect to the FUV and optical bands shows 
significant variation, being roughly 5 times greater for the 
Trapezium stars than for the exciting star of M43. We suggest 
that it is the presence or otherwise of variations in the broad 
band illuminating spectrum that is the principal determinant of 
the observed spatial variations in intensity ratios. 

4.3. The importance of the constant equivalent width 

The observed constant value of EW([Ni],Corr), coupled 
with the lack of variation in (/lL^)Fuv/(^i',))vis implies (equa- 
tion [B5]) that the ratio of [Ni] albedo to dust-scattering 
albedo is also constant within and between the nebulae: 
m,,„/m,,,, = (4 + 2) X 10"'*. The Cloudy model of § 3 im- 
plies that ^5199 ^ 2 X 10 for a non-thermal broadening of 



4.4. Variation in the /([Ni])//(Hj6) ratio 

The ratio /([Ni])//(Hy6) (Fig. la) increases by roughly a 
factor 3 between the inner and outer regions of M42 and by 
a factor of 4 to 8 between M42 and M43. The different val- 
ues of (/lL^)puv/('li'/i)Euv for the stellar populations (Table 2) 
can fully explain the difference between M42 and M43 but 
can only account for half of the variation within M42 and 
cannot explain any of the variation within M43 (where the 
single dominant star means that the illuminating spectrum 
should be constant). This implies that a small systematic in- 
crease with radius within each nebula of the ratio of albedos 
Wsiw/wh/s may also play a role. The analysis of § B.2. 1 shows 
that wiys ^0.1 when the illumination and viewing angles are 
face-on, which, combined with the above value of mstm and 
using equation (B3), implies /([Ni])//(H/3) ^ 0.03 for illu- 
mination by the Trapezium spectrum, as is observed for the 
innermost regions of M42. Inspection of Figure 1 1 shows that, 
as long as the plane-parallel approximation is valid, variations 
in the viewing angle cannot account for the inferred increase in 
Wsito/wHjS with radius since the two albedos depend on angle 
in a similar way, except for close to edge-on orientations where 
Wsito/wHjS is predicted to decrease. However, as discussed in 
§ B.4, a finite curvature of the scattering layer has the effect 
of limiting the limb brightening for edge-on viewing angles, 
and this effect is much greater for Hfi, where the scattering 
layer is much thicker than for [Ni]. This effect may explain 
the increase in vJsm/vjHfi if the average viewing angle became 
increasingly edge-on in the outskirts of the nebula (see § B.3 
for further discussion). An alternative explanation could be 
an increase with radius of the turbulent broadening within the 
fluorescent [N i] layer, but there is no independent evidence for 
such an increase. 

4.5. The /([N i])//(Hy6) ratio in other classes of nebulae 

The original motivation for this work was to calibrate mod- 
els of the formation of [N i] lines in the relatively quiescent 
Orion environment, as a step towards understanding what these 
lines indicate in the more exotic environments where they are 
unusually strong. Although continuum fluorescent excitation 
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does account for the [N i] lines in Orion, the process cannot 
produce the stronger [N i] emission seen in planetary nebulae 
or the Crab supernova remnant. 

The last two rows of Table 2 show the continuum inten- 
sity ratios produced by the SEDs of the Ring Nebula (a 
T = 1.2 X 10^ K Ranch stellai" atmosphere, O'Dell et al. 
2007) and the Crab Nebula (Davidson & Fesen 1985). The 
[Ni]/Hy6 intensity ratio scales with the FUV/EUV contin- 
uum ratio. Table 2 shows that this ratio is 7 and 3 times 
smaller for the Ring and Crab Nebulae than in Orion. Accord- 
ingly the continuum fluorescent excitation contribution to the 
[N i]/Hy6 intensity ratio produced by fluorescence will be of 
order /([Ni])//(H/3) ~ 10"-'. The observed line ratio is several 
orders of magnitude larger, showing that other processes must 
be at work. Thermal excitation by warm gas, perhaps produced 
by penetrating energetic photons or ionizing particles, seems 
to be needed. High resolution observations of the lines given 
in Figure 9 could test whether thermal processes do account 
for the observed spectrum. 

5. CONCLUSIONS 

There are multiple conclusions that can be drawn from the 
work reported upon in this paper. Some are positive, in the 
sense that we find quantitative explanations for detailed obser- 
vations of the Orion Nebula, while some are negative, in the 
sense that we demonstrate that the FUV pumping mechanism 
cannot be the dominant process in objects like the planetary 
nebula the Ring Nebula and supernova remnants like the Crab 
Nebula. The specific conclusions are summarized below. 

1 . The [N i] doublet is produced not by collisional excita- 
tion out of the lower-lying ground-state of neutral nitro- 
gen. Rather, it is the result of FUV continuum radiation 
being absorbed and populating a higher electronic state 
which then populates the upper states of the [N i] doublet 
by cascade. 

2. The process operates in the thin transition boundary of 
the PDR that is close to the overlying ionization front. 

3. This process means that one cannot use the relative 
strength of the two members of the [Ni] doublet as 
density indicators. 

4. In order for this mechanism to produce the intensity of 
the [N i] emission seen in the Huygens Region of the 
Orion Nebula there must be a non-thermal component 
to the broadening of the FUV absorption line that drives 
the process, with FWHM of approximately 5 km s ' 
(see Fig. 5). We argue in Appendix C that the origin of 
this broadening cannot be the same as the transonic tur- 
bulence that is believed to be responsible for broadening 
the optical emission lines in the H ii region because the 
latter operates at too large a scale to affect the radiative 
transfer in the thin pumping layer. Instead, we suggest 
that small-scale thermal instabilities may be responsible. 

5. The constant value of the equivalent width of the [N i] 
doublet with respect to the underlying scattered light 
continuum can be interpreted as the PDR being opti- 
cally thick to scattered starlight and a combination of 
reasonable assumptions about the scattering properties 
of the solid particles in the PDR and the orientation of 
the PDR. 



6. The efficacy of this pumping process is critically de- 
pendent upon the ratio of FUV/EUV radiation from the 
illuminating sources. We show that the stars associated 
with the Orion Nebula and the independent low ioniza- 
tion H II region M43 explain the different amounts of 
[N i] excess emission in these very different objects. 

7. The FUV/EUV ratio for a bright planetary nebula (the 
Ring Nebula) and the well observed Crab Nebula super- 
nova remnant indicate that the FUV pumping mecha- 
nism that explains the Orion Nebula and M 43 is not the 
source of the excess [Ni] emission in those objects. 
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APPENDIX 



A. THE N I EMISSION MODEL 

Here we describe recent improvements in the treatment of 
N 1 emission in the spectral simulation code Cloudy. Our model 
includes many emission processes because it is intended to be 
general, and applicable to other environments. 

A. 1 . The atomic model 

In order to optimize the speed of the model, we have chosen 
to model the N i atom using a 5 -level atom for the metastable 
levels. The fluorescence processes discussed in this paper 
(as well as the recombination pumping) are added as rates 
populating the various excited metastable levels. The level 
energies were obtained from Moore (1975) and the lowest 
five levels are listed in Table 3 and shown in Figure 7. An 
additional ten FUV lines can absorb photons in the 95 lA- 
1 161 A range. These lines drive the fluorescence process and 
will be discussed in more detail below. 

Transition probabilities have been computed by Butler & 
Zeippen (1984), Godefroid & Fischer (1984), Hibbert et al. 
(1991), Tachiev Sl Froese Fischer (2002) and Froese Fischer Sl 
Tachiev (2004). 



Table 3 

[Nl] Energy Levels. 



Configuration 


Term 


J 


Energy (cm ' ) 






3/2 


0.000 




^D" 


5/2 


19 224.464 






3/2 


19 233.177 




2 po 


1/2 


28 838.920 






3/2 


28 839.306 



35000 



30000 



25000 



0) 

E 

c 

> 
re 



20000 



5000 



U466.49 



X3466.54 



X10398 
XI 0397 

X10407 




3/2 
5/2 



X5198 
X,5200 



Figure 7. The lowest five levels of the N I model. The 5 198, 5200A pair of 
lines are denoted by /i5199^ in the text. A fourth line of the IR multiplet, 
^Pl/2 - ^D3/2 at 10408A, is not shown for clarity. 



Table 4 compares the Butler & Zeippen (1984) and Gode- 
froid & Fischer (1984) rates for the forbidden transitions, re- 
ferred to as the "1984" rates, with the more recent calculation 
of Froese Fischer & Tachiev (2004), referred to as the "2004" 
rates. The latter rates are used. Hibbert et al. (1991) do not 
give transition probabilities for the forbidden transitions. 

The electron collision rates for [N i] have been the subject 
of a number of studies. Berrington et al. (1975) computed 
electron collision cross sections which Dopita et al. (1976) 
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Table 4 

[N I] transition probabilities 



air wavelength 


transition 


1984 


2004 


5200.3 


2 no ^ 
^5/2 




ro 
■^3/2 


5.77(-6) 


7.57(-6) 


5197.9 


2 no _ 
^3/2 




ro 
■^3/2 


2.26(-5) 


2.03(-5) 


3466.543 


2 po 
'^1/2 




ro 
■^3/2 


2.52(-3) 


2.61(-3) 


3466.497 


2 po 

3/2 


^4 


ro 
■^3/2 


6.21 (-3) 


6.50(-3) 


10398.2 


2 po 

1/2 


^2 


^^5/2 


3.03(-2) 


3.45(-2) 


10397.7 


2 po 
^3/2 


^2 


^^5/2 


5.39(-2) 


6.14(-2) 


10407.6 


2 po 

1/2 


^2 


^^3/2 


4.63(-2) 


5.27(-2) 


10407.2 


2 po 

3/2 


_^2 


^3/2 


2.44(-2) 


2.75(-2) 



converted into collision strengths. These were later summa- 
rized by Berrington &. Burke (1981). Dopita et al. (1976) found 
some discrepancies with existing observations and speculated 
that the disagreement was due to uncertainties in the collision 
strengths. Tayal (2000) presented close-coupling calculation 
of the effective collision strengths while Tayal (2006) redid the 
calculation with significantly different results. 

Table 5 gives the history of these electron collision strengths. 
The *S° - D° collision strength affects the intensity of the col- 
lisionally excited contribution to the [N i] /15199^ lines, while 
the ^£>3^2 ~^ ^1/2 collision strength affects the density diagnos- 
tic. The Berrington & Burke (1981) and Tayal (2006) results 
are in reasonable agreement suggesting that the theoretical 
calculations have converged onto a stable value. 



Table 5 

History of [Nl] collision strengths at lO'^K 



Reference 


45.0 _2 QO 


2 no _2 no 
^3/2 -^1/2 


Berrington & Burke 1981 


0.48 


0.27 


Tayal 2000 


0.044 


3.24 


Tayal 2006 


0.561 


0.257 



We know of no rates for collisions with hydrogen atoms. 
This should be included since we expect that [N i] may form 
in shallow regions of the PDR, where n(H'^) >> n^,. Cloudy 
does include a general correction for H" collisions based on 
the He rate. The effective electron density, with this correction, 
is «j -H 1.7 X 10 (h") based on the discussion by Drawin 
(1969). 

Cloudy includes many other line formation processes in ad- 
dition to thermal coUisional impact excitation (Ferland Sl Rees 
1988). Continuum pumping of the FUV lines around lOOOA 
will be very important. We treat this as described in Ferland 
(1992) and Shaw et al. (2005). Fluorescent excitation is mainly 
produced by stellar FUV photons. Pumping will be efficient 
until the N i lines become optically thick At tiiis point they will 
have absorbed stellar photons over the Doppler width of the 
line. Below we explore how the pumping efficiency depends 
on the turbulent contributor to the line width. Other opacity 
sources will affect the strength of the pumped contributor to 
N I by removing FUV photons before they are absorbed by 
N I. The two most important opacity sources are extinction 
of the stellar radiation field by grains within the H 11 region 
and PDR, and shielding by the forest of overlapping H2 lines 
in deeper parts of the PDR. These processes are all included 
self-consistently in our calculations. 



A. 2. The fluorescence mechanism 

In strict LS coupling, transitions that change the total spin 
of the atom (called intercombination transitions) are forbidden 
and FUV pumping out of the quartet ground term could not 
eventually populate the doublet excited terms that produce the 
observed lines. However, deviations from strict LS coupling 
make it possible that a significant fraction of excitations by 
FUV photons will eventually populate the excited doublets. 
There are two possible routes: either direct excitation by an 
intercombination line from the ground state or indirect excita- 
tion by a resonance line followed by de-excitation through an 
intercombination line. In the case of N i both routes contribute. 
A complete list of the driving lines can be found in Table 6. 
For each of the driving lines we calculate a 2-level atom giving 
us the excitation rate for each of these transitions. We also 
calculated branching ratios for the cascade down from each of 
the upper levels. In these calculations we exclude the transition 
straight back to the ground state as this doesn't destroy the 
photon. Instead it can be absorbed over and over again until 
finally a different cascade from the upper level occurs (this 
neglects background opacities which will be discussed further 
down). Intercombination lines from the doublet system back to 
the quartet system are included in the cascade, but not tracked 
any further after that. This implies that routes quartet — > dou- 
blet — > quartet — > doublet are not included in the pumping 
rates. We expect the error introduced by this approximation 
to be negligible. The branching ratios were calculated using 
transition probabilities from Froese Fischer Sl Tachiev (2004). 
By combining all different routes in the cascade we could cal- 
culate a probability that an excitation of a given driving line 
would result in populating any of the metastable levels. The 
results of these calculations are shown in Table 7. The list 
of intercombination lines populating the doublet metastable 
levels after an excitation in the Indl driving line is given in 
Table 8. 

In the previous discussion we mentioned that transitions 
in any of the driving lines straight back to the ground level 
were not counted because these photons would simply be re- 
absorbed until a different cascade occurs. This assumption is 
not entirely correct as there is a finite probability P^est that the 
photon is destroyed before it can be absorbed again (e.g. due to 
background opacities such as the grain opacity or bound-free 
opacity of elements with sufficiently low ionization potentials). 
Additionally there is a probability Pgsc that the photon escapes 
from the cloud before it can be absorbed again. In order to 
account for these processes, we modify the excitation rate 72 
in s"' obtained from the 2-level atom as follows: 



1 -yS(l -^dest-^esc) 



(Al) 



where y6 is a constant that gives the fraction of excitations in a 
driving line that is followed directly by a de-excitation back to 
the ground level. ¥01 Indl {3 = 0.7955 and forD;>7 /? = 0.1384. 
For all other driving lines fi < 0.01 and is assumed to be zero. 
Given this formula we can the write the total pump rate for 
each of the metastable levels as 



7meta 



=z 



PL 



(A2) 



where the summation runs over all the driving lines and the 
constants Pp^^^p are given in Table 7 for each of the metastable 
levels and driving fines. 
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Table 6 

The lines driving the Nl fluorescence. For each of the lines, the lower level is the ground state of Nl. The level energies are taken from Moore (1975) and the 
transition probabilities from Froese Fischer & Tachiev (2004). The effective collision strength is given by the fitting formula T = exp(a + b min[ln T, 10.82]) where 

the original data were obtained from Tayal (2006). 



label 


upper level 


Ei (cm"') 


At,(s-l) 


a 


h 


Indl 


2i22p2(3P)3d 4pj^2 


104 825.110 


1.62(8) 


-11.3423 


0.8379 


Dirl 


2i22p2(3P)3c( 2pj^2 


104 810.360 


1.95(7) 


-12.3982 


0.7458 


Dir2 


105 143.710 


8.29(5) 


-9.4523 


0.3865 


Dir3 


2i22p2(3p)3d 2p3^2 


104 615.470 


4.29(5) 


-12.5580 


0.7330 


Dir4 


2i22p2(3P)4i 2p3/2 

2.«22p2(3P)3d 2p,/2 


104 221.630 


3.75(5) 


-10.8813 


0.6853 


Dir5 


104 654.030 


2.63(5) 


-13.6532 


0.7712 


Dir6 


2/2p2(3P)3fi(2i5j^2 


105 119.880 


1.71(5) 


-9.9035 


0.3919 


Dir7 


2i22p2(3P)4.«2Pi/2 


104 144.820 


1.69(5) 


-11.4470 


0.6734 


Dir8 


2i22p2(3P)3i 2p3/2 


86 220.510 


4.94(4) 


-5.4776 


0.1789 


Dir9 


2i22p2(3P)3i2p,/2 


86 137.350 


2.72(4) 


-6.3304 


0.1966 



Table 7 

Probability Ppump of populating a metastable level after an excitation in each of the driving lines. 

level Indl Dirl Dir2 Dir3 Dir4 Dir5 Dir6 Dir7 Dir8 Dir9 

2d°^2 0.0417 0.0468 0.3408 0.2328 0.7937 0.1338 0.0623 0.0238 0.6615 0.0000 

2d°^2 ''■3441 0.8621 0.0233 0.0895 0.1068 0.1644 0.2908 0.8397 0.0694 0.7369 

2p°^2 0.0113 0.0239 0.0090 0.1617 0.0167 0.4404 0.4881 0.0876 0.0450 0.1777 

2p°/2 0-0112 0.0265 0.6253 0.5108 0.0824 0.2588 0.1569 0.0484 0.2240 0.0854 



Table 8 

The list of intercombination Unes that can occur after an excitation in the Indl 
driving line. For each line the upper level is 2s^2p-{-iV)3d ''P5/2. The level 
energies are taken from Moore (1975) and the transition probabilities from 
Froese Fischer & Tachiev (2004). 



lower level E, (cm ' ) A^, (s ' ) 



2s^2p^ ^-Dili 


19 233.177 


1.24(7) 


2i22p2(3P)3p 2D3/2 


96 787.680 


2.92(6) 


2^22^3 2D5/2 
2i22p2(3P)3p 2D5/2 


19 224.464 


1.30(6) 


96 864.050 


2.16(5) 


2s^2p^ ^Viij 
2.s22p2(3P)3p 2p3/, 


28 839.306 


1.06(5) 


97 805.840 


9.87(3) 



For completeness we should also mention that pumping 
of the metastable states through recombination from is 
also included in our modeling. We use the formulas given 
in Pequignot et al. (1991). These only give the rates to the 
full and metastable terms. In our modeling we split up 
these rates for each level according to statistical weight. This 
pumping mechanism will of course only be effective inside 
the ionized region as nitrogen has a slightly higher ionization 
potential than hydrogen. 

From the data in Table 6 it is clear that all driving lines have 
wavelengths longward of the Lyman limit. This implies that 
the fluorescence mechanism is effective beyond the ionization 
front in the FDR. Since the temperature in the PDR is generally 
too low to coUisionally excite the metastable doublet states, 
fluorescence can even become the dominant excitation mecha- 
nism for the forbidden N i lines in the PDR. It should also be 
noted that even very weak direct excitation lines can have a 
significant contribution to the fluorescence mechanism. If the 
PDR has sufficient column density, then all driving photons 
will eventually be absorbed. A low transition probability in the 
driving line only means that the effect is spread over a larger 
area. 



The fluorescence mechanism will produce permitted Ni 
emission lines that are observable in deep spectra. The cascade 
routes that populate the metastable levels will produce lines in 
the doublet system with wavelengths ranging between 8567 A 
and 5.382 //m, as well as UV lines that cannot be observed 
from the ground. The shortest wavelength lines will tend to 
be the strongest since they come from the lowest levels where 
there are only a few alternative routes the cascade can take. 
In the ionized region these lines can also be produced by re- 
combination from — » N", but in the PDR these lines can 
only be produced by the fluorescence mechanism described 
here. So if these lines are observed in the PDR, it is conclusive 
proof for continuum pumping of the [Ni] lines. In Table 9 
we list all optical cascade lines with wavelengths shorter than 
1 yum and in Table 10 we list the branching probabilities for 
each of the driving lines. Excitations by the Dir8 and Dir9 
driving lines only produce UV cascade lines and are therefore 
not included in Table 10. Each driving line has its own char- 
acteristic spectrum. The relative strength of the contribution 
for each driving line depends on the incident spectrum, the 
optical depth in each of the driving lines and the escape and 
destruction probability for the Indl and Dirl driving lines. So 
no generic prediction for the spectrum can be made. However, 
a straight average indicates that the AA 9387, 9029, 9060, and 
8629 lines will be the strongest, with the A 9387 line having 
about 1.8% of the flux of the /15199+ doublet. It should be 
noted that the cascade lines shown in Table 9 are not predicted 
by Cloudy. 

At this point we should discuss the accuracy of the transition 
probabilities of the intercombination lines. Accurate values for 
such lines are hard to obtain since they are quite sensitive to the 
details of the calculation. However, for direct excitation lines, 
accurate values for the transition probability are not crucial. 
Using the argument from the previous paragraph it becomes 
clear that an error in the transition probability would only 
imply that the absorption of the driving photons would happen 
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Table 9 

The list of permitted lines in the doublet system that can be excited by the 

fluorescence mechanism described here. Only lines with wavelengths 
between 8567 A and 1 /jm are listed. The transition probabilities were taken 
from Froese Fischer & Tachiev (2004). 





transition 






A*,(s-') 


OJO / . IJJ 


9 (.29 „2/op\'i „ 


2po 

'^3/2 


i 9 .29 rt^j'^PA'i c. 2p, _ 


A RUft^ 
^■.o /yO) 


8594.000 


2s^2pH3P)3p 


2po 


-> 2i22p2(3P)3i2pi/2 


2.10(7) 


8629.235 


2s^2p\3P)3p 


2po 


^ 2i22p2(3P)3.v 2p3/2 


2.68(7) 


8655.878 


2s^2pH3P)3p 


2po 


^ 2i22p2(3P)3l ^P3/2 


1.08(7) 


9028.922 


2s^2p^{3P)3d 


^Pl/2 


2s^2p^(3P)3p -S°/2 


3.20(7) 


9060.475 


2s^2p^i3P)3d 


^P3/2 


^ 2s^2p^(3P)3p 2S°/2 


3.21(7) 


9386.805 


2s^2p\3P)3p 


^3/2 


-» 2/2p2(3P)3i2p|/2 


2.14(7) 


9392.793 


2s^2pH3P)3p 


'D5/2 


^ 2i22p2(3P)3i ■^-Pin 


2.52(7) 


9395.848 


2s^2p^{3P)4s 


^P3/2 


2.v22p2(3P)3p 2s°/2 


1.81(4) 


9460.676 


2s^2pH3P)3p 


^3/2 


^ 2s^2p^{3P)3s 2P3/, 


3.74(6) 


9464.169 


2s^2p^{3P)4s 


^Pl/2 


2i22p2(3p)3p 2s°/2 


3.50(5) 



over a smaller or larger area, but the total amount of pumping 
would remain the same when integrated over the entire PDR. 
This of course assumes that the PDR is optically thick. If 
that is not the case, then an error in the transition probability 
would alter the escape probability of the driving line. In such 
circumstances accurate transition probabilities are needed. For 
indirect excitation lines, accurate transition probabilities are 
always needed (even when the PDR is optically thick) since 
the intercombination line has to compete with stronger, fully 
allowed transitions in the cascade down from the upper level of 
the driving line. However, this problem is mitigated by the fact 
that there is only one indirect driving line versus nine direct 
driving lines. So an error in this component would only have a 
limited effect on the total pumping. In Table 1 1 we compare 
the transition probabilities of the lines involved in the cascade 
down from the indirect excitation using data from Hibbert et al. 
(1991, length form) and Froese Fischer & Tachiev (2004). It 
is apparent that discrepancies up to 1 dex and more can occur, 
indicating the difficulty in calculating these data. 

A. 3. Density - temperature diagnostics in the collisional 
excitation case 

If the lines were collisionally excited then the electron den- 
sity could be determined from the ratios of the intensities of 
two lines of the same ion, emitted by different levels with 
nearly the same excitation energy (AGN3). Temperature is 
indicated by emission from levels with different excitation en- 
ergies. Together, the gas pressure could be directly measured. 
This can test whether the lines are thermally excited, and is 
useful for reference by future studies which will look into the 
formation of [N i] lines in planetary nebulae, the Crab Nebula, 
and cool-core cluster filaments. 

We show several emission line diagnostics for the collision- 
ally dominated case, using our updated atomic data and model 
atom. Line pairs such as the ratio of lines of [N i], 

R„ = I{^Dsi2 ^S3,2)ll(^Di,2 ^ ^83/2) = .15200/^5198 

(A3) 

indicate the electron density in gaseous nebulae, as shown by 
Seaton & Osterbrock (1957) and Saraph & Seaton (1970) for 
[0 11]. Note that the energy order of the J levels within the 
term depends on both the charge and electronic configuration. 
The line ratio is defined so that it decreases as density increases. 



Every collisional excitation is followed by the emission of a 
photon in the low-density limit. Since the relative excitation 
rates of the ^^5/2 and ^1)3/2 levels are proportional to their 
collision strengths, the ratio is 



R„(ng — > 0) 



reosn -^83/2) ^ 0.337 

T(2Z)3/2-4 53/2) 0.224 



1.5 



(A4) 



This is valid when kT » 6e, where Se is the difference in ener- 
gies of the upper levels. This holds for all temperatures where 
the optical lines emit due to the small energy difference of the 
upper levels. In the high-density limit collisional processes 
dominate and set up a Boltzmann level population distribution. 
The relative populations of the ^Ds/2 and ^0^/2 levels are in 
the ratio of their statistical weights, and the relative intensities 
of the two lines are in the ratio 



^ ^^('£'5/2)A,5200 3 7.57 x 10"^ ^ 

R„ («,. — > 00) = — = — 0.60 

t^(2D3/2)A,5198 2 2.03 X 10-5 

(A5) 

The line ratio varies between these intensity limits as the den- 
sity varies. The critical density, the density where the col- 
lisional and radiative deexcitation rates are equal, is ncWr ~ 
10^ cm-3 at ~ 10"* K. 

Figure 8 compares the [N i] R„ with the more commonly 
used [On], [S 11], and [Clm] density indicators using data from 
B2000. The [0 11] collision strengths computed by Kisielius 
et al. (2009) were used. The behavior of these curves is quali- 
tatively similar, going to the ratio of statistical weights at low 
densities and a ratio that depends on the radiative transition 
probabilities at high densities. 



00 



I 




0.3 



log/7. 

Figure 8. [N l] and three of the commonly-used density indicators present in 
optical spectra. The points are from B2000. 

The gas kinetic temperature can be determined from ratios 
of intensities from two levels with considerably different exci- 
tation energies. For [Ni] we have the ratio. 



= /(2p ^4 S)ll(^D 

(3466.49 + 3466.54)/(5198 + 5200) 



=^3467+ M5 199+ 



(A6) 

These density-temperature indicators can be combined to 
form a unified diagnostic diagram, as has long been done for 
[Oiii] (AGN3 Figure 5.12). Figure 9 shows calculated curves 
of the values of the two [N i] intensity ratios for various values 
of T and n^. 
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Table 10 

Branching probabilities of the cascade lines for each of the driving lines. Excitations by the Dir8 and Dir9 driving lines do not produce any optical cascade lines. 



J - 


Indl 


Dirl 


Dir2 


Dir3 


Dir4 


Dir5 


Dir6 


Dir7 


ojO / . / jD 


n nnnn 
u.uuuu 


n nnn i 

U.UL/U 1 




ft nm 


U.UUo J 


n on 1 s 


U.UUj J 


U.UUZJ 




(1 nono 


0000 


0000 


0.0031 


008? 


0.0129 


0.0549 


0.0229 


8629.235 


0.0002 


0.0008 


0.0868 


0.0086 


0.0469 


0.0082 


0.0190 


0.0135 


8655.878 


0.0000 


0.0000 


0.0000 


0.0016 


0.0041 


0.0065 


0.0275 


0.0115 


9028.922 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.2840 


0.0000 


0.0000 


9060.475 


0.0000 


0.0000 


0.0000 


0.2546 


0.0000 


0.0000 


0.0000 


0.0000 


9386.805 


0.0597 


0.1270 


0.0004 


0.0045 


0.0040 


0.0196 


0.0343 


0.0532 


9392.793 


0.0052 


0.0058 


0.0487 


0.0256 


0.0534 


0.0000 


0.0059 


0.0000 


9395.848 


0.0000 


0.0000 


0.0000 


0.0000 


0.0002 


0.0000 


0.0000 


0.0000 


9460.676 


0.0104 


0.0222 


0.0001 


0.0008 


0.0007 


0.0034 


0.0060 


0.0093 


9464.169 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0027 



Table 11 

Comparison of the transition probability for various intercombination lines 
used in our model. For each line the upper level is 2s^2p^(3P)3d ''P5/2. The 
transition probabilities (units s"') are taken from Hibbert et al. (1991, labeled 
"1991") and Froese Fischer & Tachiev (2004, labeled "2004"). The latter were 
used in our model. 



lower level 



1991 



2004 



2^22^3 2d3^2 
2/2p2(3P)3p 2Dvo 
2^22^3 2D5/2 
2/2p2(3P)3p 2D5/2 

2.v22p2(3P)3p ^Py2 



9.08(5) 
3.43(5) 
1.02(6) 
6.48(4) 
1.04(6) 
5.67(4) 



1.24(7) 
2.92(6) 
1.30(6) 
2.16(5) 
1.06(5) 
9.87(3) 
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Figure 9. The derived density and temperature, in units of log n^ (cm ) and 
10* K respectively, as deduced from line intensity ratios from the model [N I] 
atom. 



We used the line intensities measured in bright central re- 
gions (B2000; Esteban et al. 2004, hereafter E2004) to estimate 
Hp and T. B2000 presented high-resolution spectrophotometric 
observations of the Orion Nebula in the AA3500 - 7060A range. 
Their slit position was 37" west of 0' Ori C. This is close to 
the position modeled by Baldwin et al. (1991). E2004 covered 
the /1/13100 - 10400A range. Their slit was oriented east-west 
and centered at 15 arcsec south and 10 arcsec west of 0' Ori C. 



The value of R„ was 0.60 ± 0.02 for the average of the 
blue and red spectra in B2000 and 0.59 for E2004. This is 
plotted in both Figures 8 and 9. The results are surprising — 
the electron density indicated by the [N i] lines, which should 
form in partially ionized gas giving lower electron densities, 
is 0.2-0.4 dex larger than densities indicated by lines which 
form in highly ionized regions. 

It is not now possible to measure the temperature using 
the /13467^ line. We know of no detection of this line in 
the Orion environment. However Esteban et al. (1999) re- 
ported the upper limit of /(3467+)//(Hj0) < 10"^ from the data 
presented by Osterbrock et al. (1992). This corresponds to 
i3467+//l5199+ < 0.22. 

These values of the line ratios are shown in Figure 9. The 
temperature limit indicated by the line is consistent with forma- 
tion in a photoionized environment. The high density would 
be surprising if true. Actually this can be taken as independent 
evidence that the lines do not form by collisional excitation. 

A.4. Dissociation of nitrogen-bearing molecules 

Storzer & Hollenbach (2000) show that significant optical 
[O i] emission can result from dissociation of oxygen-bearing 
molecules. Could an analogous process contribute to the [Ni] 
emission we observe in Orion? 

Storzer & Hollenbach (2000) consider OH photodissociation 
and subsequent [O i] 6300-1- emission in detail. The intensity of 
the line that is produced depends on the photodissociation rate, 
the branching ratio for populating the excited level producing 
a particular line, and the extinction between the point where 
the emission is produced and the surface of the cloud. 

We include molecular photodissociation by the processes 
included in a modified version of the UMIST (Le Teuff et al. 
2000) database (Rollig et al. 2007). Our original treatment, 
described in Abel et al. (2005), considered each reaction on 
an ad hoc basis. In our upcoming release we will generalize 
our treatment of the chemistry to more systemically consider 
reactions, their inverses, and maintain an accounting of the 
consistency (Williams et al., in preparation). This is a step 
towards treating the chemistry as a coupled system that is 
driven by external databases. 

Using the results from the chemistry network we can identify 
all photodissociation processes. The N-bearing molecules NH, 
CN, N2, NO, and NS produce N*^ following photodissociation. 
We save this photodissociation rate per unit volume at each 
point in the cloud, and assume that each dissociation produces 
N*^ in the ^D" level. Using that we can estimate the contribution 
of this pumping process to the production of the /15199^ lines. 



Pumping up the [N i] nebular lines 



13 



The observed emission is predicted by attenuating the local 
emission by the absorption optical depth from the creation 
point to either side of the cloud. Grains are the dominant 
opacity source at optical wavelengths for conditions similar to 
the Orion Nebula. This produces an upper limit to the emission 
because of the assumption that 100% of photodissociations 
produce in the excited state producing [Ni] 5199^. This 
upper limit is added to the flux of the /15199^ lines. 

Figure 4 shows that there are regions of the cloud where 
photodissociation could make [N i] emission. However this 
is deep enough within the cloud that the process makes no 
significant contribution to the observed flux. The process may 
be important in other environments, however. 

This physics is included in the current release of Cloudy 
(ClO.OO). 

B. EFFECTIVE ALBEDOS FOR GENERALIZED 
SCATTERING PROCESSES 

In this appendix, we develop a simple framework for eval- 
uating the intensity of radiatively driven continuum and line 
processes in a photoionized nebula, which will elucidate the 
dependence of line ratios and equivalent widths on the shape of 
the exciting stellar spectrum and on geometrc factors. All three 
emission mechanisms, dust continuum, [N i], and Hj6 lines, can 
be thought of as diflFuse reflection or scattering processes in 
the broadest sense, with each being driven by a diflferent wave- 
length band of the stellar continuum:^ 

• Dust continuum is coherent scattering in the optical 
sense of visual band photons (~ 5000 A). 

• Fluorescent [N i] is highly incoherent scattering of FUV 
pumping photons (~ 1000 A) into visual photons. 

• To the degree that static photoionization equilibrium 
holds, then HjS emission is scattering (albeit in a sta- 
tistical and indirect way) of ionizing EUV photons 
(< 912 A) into visual photons. 

In each case, one can define an effective albedo m, which is an 
efficiency factor that relates the intensity of scattered or emitted 
photons to the intensity of incident photons (see Figure 10 and 
§ B.l below). Therefore, any variation in the observed intensity 
ratios must be due to either (1) variations in the spectral energy 
distribution (SED) of the stellar radiation field, or (2) variation 
in the effective albedos, or (3) a breakdown of the simplifying 
assumption of a single infinite plane-parallel scattering layer. 
In the remainder of this appendix we discuss in detail the 
contributions of (2) and (3), while the role of (1) is explored in 
§4.1 above. 

B.l. Formal calculation of intensity ratios and equivalent 

widths 

Under this black-box "scattering" or "reprocessing" descrip- 
tion, the efficiency of the scattering can be described by an 
effective albedo m^s (see Fig. 10), so that the photon intensity 
/ for each line or continuum process is proportional to the 
local continuum flux Fq in the spectral band that excites the 
scattering process: 



pQ^eS , , -1-2 -1 

photons s cm sr 



(Bl) 



' The [O I] 6300 A emission is more complicated because of a strong 
dependence on the ionization parameter, and so will not be considered further 
here. 



Scattering layer 




Incident radiation 



Figure 10. A black-box approach to generalized scattering processes. Mono- 
directional radiation with a flux parallel to its beam Fo (photons s"' cm"^) 
is incident on a plane parallel scattering layer from a direction = cos Bg, 
where Bq is the angle from the normal to the layer. The azimuth of the incident 
radiation may be taken as 0o = without loss of generality. The intensity of 
emergent scattered radiation in a direction /i, ^6 is / (photons s"' cm"^ sr"' ), 
where the scattered radiation may be in a very different wavelength band 
from the incident radiation. The effective albedo of the scattering process will 
depend on the directions of both the incident and emergent radiation and is 
defined as m(i.io; fi, </>) = AttI/Fq. 



where 



^ {'^i^)band(A'^)band , , -1 -2 

Fi) = r photons s cm . 



(B2) 



In this expression, R is the distance from the star to the scat- 
tering layer, and (/lL^)band and (A/l)band are respectively the 
mean SED and wavelength width of the continuum band that 
excites the process. In general, the albedo will be a function 
of the illumination angle and viewing angle (Fig. 10). The 
particular values of the albedo for the production of the H/3 re- 
combination line, the [N i] fluorescent lines, and dust-scattered 
continuum in the nebula are calculated in § B.2 below. 

The line ratios and equivalent widths measured in § 2 will 
then be given by 



/([Ni]) ti7_„„<^L,)p,v(Ai) F 



E 



EW(I^, Corr) = 



EW([Ni],Con-) 



mufi </1L^)euv (A/1)eu 

TTTsm (/tL^i)Fuv (A/{)fi 



(B3) 



(B4) 



(B5) 



A particularly simple limiting case is provided by the situ- 
ation in the extreme outskirts of the Extended Orion Nebula. 
For these regions studies have shown that, except for the low- 
est ionization lines, essentially all the radiation, including the 
emission lines, is scattered by dust rather than being emitted 
locally (O'Dell & Goss 2009; O'Dell & Harris 2010). For 
positions far outside the bright core of the nebula, it is rea- 
sonable to make the additional assumption that the angular 
distribution of the incident radiation (as seen by the scatterers) 
is on average similar for the continuum (which comes from 
the star cluster) and for the emission lines (which come from 
the nebular gas). That being the case, all geometrical factors 
will cancel out and the effective albedo for an emission line 
will be the same as that for the adjacent continuum, so that the 
equivalent width of the line will be simply EW(Int) - L\ine/L^, 
where Lune is the total intrinsic line luminosity of the nebula 
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and L,i is the total intrinsic continuum luminosity of the star 
cluster.^ One would therefore expect that the observed cor- 
rected equivalent widths should tend to a constant value of 
EW(Int) in the extreme outskirts of the nebula. Just such a 
behavior is seen in the observed H/? equivalent width (O'Dell 
& Harris 2010, Fig. 8), which around the outer rims of the 
nebulae tends to a value of ~ 150 A for M42 and ~ 100 A 
for M43. Comparison with the intrinsic H/3 equivalent widths 
in Table 2 shows good agreement in the case of the M43, al- 
though for M42 the predicted value of 213 A is rather higher 
than is observed. 

In Table 2 we show SED ratios between different wavebands 
for OB stellar populations characteristic of the inner Orion 
nebula (Trapezium region), the outer Orion nebula, and M43. 
These can be compared with different observed emission ra- 
tios shown in Figures 1 and 6: FUV/visual corresponds to 
EW([Ni],Coit); EUV/visual corresponds to EW(Hj0,Corr); 
FUV/EUV corresponds to /([Ni])//(H^).'^ 

B.2. Estimation of effective albedos for particular processes 
B.2.1. Ufi recombination line 

Assuming a thin, plane-parallel, ionization-bounded layer of 
dust-free hydrogen that is illuminated by ionizing photons with 
a flux Feuv (photons s ' cm"^) incident from a direction yUo, 
then the condition of global static photoionization equilibrium 
is given by 



-I 



ffB npHe dz, 



(B6) 



in which q-b is the "Case B" recombination coefficient and «p, 
He are the proton and electron densities. At the same time, the 



emergent intensity / (photons s cm sr ) of the recombi 
nation line Hfi is given by 



=- r 



oh/s n^ris dz. 



(B7) 



where a-ap is an effective recombination coefficient that only 
includes those recombinations that give rise to the emission 
of an HyS photon. Combining these, the effective albedo (see 
Fig. 10) is found to be 



Q'H/3 \ yUO 



(B8) 



where {aiipla^ ^0.12 for typical H n region conditions (Os- 
terbrock & Ferland 2006). 

The correction to this result for the presence of helium will 
be small, but the presence of dust in the ionized gas may have 
a much larger effect. Both the incident ionizing radiation and 
the emergent emission line will be affected by dust absorp- 
tion. The fraction of ionizing photons that are absorbed 
by dust is an increasing fraction of the ionization parameter 
(Aannestad 1989; Arthur et al. 2004), but for the conditions 
found in Orion reaches a maximum value of about 30% so long 
as the illumination is close to face-on. The effect is greater 
for edge-on illumination, but such cases, with ^ 1. have 
already a small albedo and so will contribute little to the ob- 
served emission so long as a variety of illumination angles is 

' As described in § 2, the contribution to the observed equivalent widths of 
the atomic continuum emission from the nebular gas should first be corrected 
for. 

* Note that only 2 of these 3 quantities are independent. 



present (see discussion in § B.3). The absorption of emergent 
Mfi photons is more important since this is largest for precisely 
those cases yu «; 1 which would give the highest albedo in the 
dust-free case (eq. [B8]). If the dust absorption optical depth of 
the scattering layer at the wavelength of H/J is r, then in the ap- 
proximation that the ionized density is constant equation (B8) 
becomes 



CTHyS = (1 



\QfB/ T ^ 



(B9) 



The maximum relative boost in the albedo due to limb bright- 
ening as — > 0, which is infinite in the dust-free case, is now 
limited to 1/ (1 - e^'^), which is a factor of 3-5 for the values 
of T ^ 0.1-0.3 expected in Orion. Note that scattering by dust 
of the Hy6 photons is ignored in this approximation. 

B.2.2. Fluorescent 

The FUV fluorescent pumping of the optical [N i] lines is 
only efficient at wavelengths where the opacity of the pump- 
ing line exceeds the background continuum opacity, which at 
FUV wavelengths is dominated by dust. This gives a limit 
6A — ASvjc to the wavelength interval that contributes to the 
pumping, where 6v is of order the Doppler width of the line.^ 
If there are a number A^une pumping lines, each of effective 
width 5v, then the fraction of the total FUV continuum that 
contributes to the pumping is ^ Mine(5t;/c)((/l)puv/(A/l)Fuv), 
where (/1)fuv is the average wavelength of the pumping lines 
and (A/l)puv is the wavelength width of the FUV band.'" If a 
fraction f„g of all pumps results in the emission of a line in the 
optical [Ni] /1/15 198, 5200 A doublet, then the effective albedo 
for "scattering" of FUV continuum into these lines is 



47r/„ 



~ 75199-^line j 



6v\l {A), 



(A/l)p, 



-T"Fuv//^o ^-Tsm/f^ 

(BIO) 



where Tpuv and r,,,, are the continuum absorption optical depths 
between the star and the pumping layer, measured perpendicu- 
lar to the layer, and at the wavelengths of the pumping FUV 
lines and the emerging optical lines, respectively. 

The opacity of the FUV pumping lines rpump is proportional 
to the abundance of N", and so is very low inside the ionized 
gas, rising suddenly at the ionization front. Therefore, for 
strong pumping lines, the fluorescent excitation (which peaks 
at Tpump - l-io), is concentrated in a thin layer just behind the 
ionization front so that and r^m are insensitive to variations 
in the illumination cosine fiQ. For the weakest pumping lines, 
on the other hand, the pumping layer extends deeper into the 
neutral PDR and so Tpuv and Tj,,, are generally larger and 
become roughly proportional to yUo. 

Fi gure lib shows results for zusm for the case of per- 
pendicular illumination /iq = 1 and assuming Mine = 10, 
6v = lOkms-i, /„„ = 0.1, </1>puv/(A/1)fuv = 4.3, and 
Tpuv = 1.5t5,5,. It can be seen that the same optical depth of 
dust has a considerably larger effect on the [N i] albedo than on 
the HjS albedo, particularly for oblique viewing angles (small 

' For a Gaussian line profile, dv = b Vln(^o/^dusi), where b = 0.601 X 
FWHM is the Doppler broadening parameter, is opacity at line center, and 
is the continuum dust opacity. 
For the wavelength range of 950-1200 A used in Table 2, 
<-1>fuv/(A^)fuv = 4.3. 
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Figure 11. (ci) Effective albedo (cth,?, Eq. [B9]) for "scattering" by a plane ionized layer of normally incident ionizing EUV photons into optical HyS photons as 
a function of the viewing direction /j. Results are shown for three different values of r the perpendicular dust absorption optical depth through the layer at the 
wavelength of liJ3. The dust acts primarily to limit the limb brightening at small values of /i. Note that mnp is defined in terms of numbers of photons; in terms of 
energy ttie values would be a factor /ivh/j/(/!v>euv -0.15 times smaller, {h) Same as (a), but for fluorescent "scattering" of incident FUV photons into optical [N I] 
photons {rjTsm, Eq. [BIO]). 



/j). This is because the dust absorption layer completely over- 
lies the fluorescent scattering layer in the [N i] case, whereas 
in the case of Hfi the dust is mixed in with the line-emitting 
gas. As a result, whereas the H/3 albedo wu/i simply saturates 
at small the [N i] albedo ZtJsigg has a maximum at p — ^5199 
and then drops to zero as yU — > 0. 

B . 2 . 3 . Scattered starlight 

The dust scattering of starlight in the nebula can be divided 
into two parts: (1) back-scattering by dust in the PDR and 
molecular cloud located behind the the nebula, which has a 
high optical depth, and (2) small-angle scattering by dust lo- 
cated in the diffiise clouds in front of the nebula (the neutral 
veil, Abel et al. 2004, 2006), which has a smaller optical depth 
(r = 0.1-1, O'Dell & Yusef-Zadeh 2000). In both cases, the 
results will be sensitive to the optical properties of the dust 
grains, which at the simplest level can be characterized by the 
single-scattering albedo tuq, which is the probability that a pho- 
ton interacting with a grain is scattered rather than absorbed, 
and the asymmetry parameter g, which is the mean cosine of 
the scattering angle {g - for isotropic scattering). Dust in 
Orion is found to have a high value of the total/selective extinc- 
tion ratio Ry - 5, possibly due to grain coagulation (Cardelli 
& Clayton 1988). Theoretical calculations of the optical prop- 
erties of a grain population with this value of Ry (Fig. 4 of 
Draine 2003) imply that at optical wavelengths (~ 5000 A) 
the albedo is relatively high {■utq ^ 0.8) and the scattering 
is moderately forward-throwing {g ^ 0.6), whereas at FUV 
wavelengths (~ 1000 A) the albedo is lower (nio - 0.4) and 
the scattering is extremely forward-throwing (g ^ 0.8). Obser- 
vations in Orion of scattered FUV continuum (Shalima et al. 
2006) and scattered optical emission lines (§ 3.1 of Henney 
1998) are consistent with these values, although in both cases 
it is only a combination of -utq and g that is constrained. Earlier 
studies (Schiff'er & Mathis 1974; Mathis et al. 1981; Patriarchi 
& Perinotto 1985) have found different values, and even evi- 
dence that the dust properties vary with position, but the results 
are very sensitive to the assumed geometry of the scattering. 
In the following we present results for both high-albedo and 
low-albedo grains, which can be taken as representative of the 
range of possible optical grain properties at visual and FUV 
wavelengths. 



The problem of back-scattering by the molecular cloud has a 
well-known solution in the case of isotropic scattering (Chan- 
drasekhar 1960), giving an effective albedo of 

T^... = H(n) Hi^io) (B 1 1) 

where H(p) is the Chandrasekhar H-function. An approximate 
analytic form for the H-function (Henney 1998), accurate to 
< 5% for mo < 0.9 is 

mp) = 1 -i-0.5cTo/i(l + 1.8;[/"'*tir^)ln(l (B12) 

For the more relevant case of asymmetric scattering {g + 0), 
the problem is more difficult to solve, and is no longer axially 
symmetric unless yUo = 1 ■ However, for illumination angles 
that are not far from face-on, a good approximation is found by 
simply multiplying the isotropic results by a factor of (1 - g)^^ 
(Henney 1998). The results of this approximation are shown 
in Figure 12fl, where it is seen that typical values of tiTj„„ — 
0.2-0.3 are obtained at visual wavelengths, but much smaller 
values (tiTj„„ < 0.05) are seen at FUV wavelengths. In both 
cases, the scattering is approximately Lambertian (brightness 
independent of viewing angle) when the illumination is close 
to face on. Note however, that this approximation ignores 
the fact that as yUo is decreased, then the forward-throwing 
part of the phase function begins to be sampled at small for 
favorable viewing azimuths 0, which would tend to increase 
the limb brightening for <\. 

For the case of forward scattering, one can use the results 
for diff'use transmission through a homogeneous plane parallel 
layer of optical thickness t (Chandrasekhar 1960), where the 
eff'ective albedo can be expressed in terms of Chandrasekhar's 
X and Y functions: 

Wj„„ = yt/o (yU - pt-iif^ Wo 0(//,yU(), ^) 

X [y(yU,T)X(A^o,T) - X(//,T)y(yUo,T)] . (B13) 

where <1) is the scattering phase function and X and Y depend 
implicitly on and tuq. In the limit of small t, it is sufficient to 
include only single scattering, which yields the approximation 
X*'*(yU, r) = 1, y*'*(//, r) = e^'^^^. For multiple scattering in 
the isotropic case, extensive tables have been published for 
X and Y (e.g., Mayers 1962) and we find that an acceptable 
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Figure 12. (a) Effective albedos for back-scattering from a very optically thick dusty layer as a function of the viewing angle The upper two lines (solid and 
dashed) are for optical properties typical of Orion dust at visual wavelengths, while the lower two lines (dotted and dot-dashed) are for those of the same dust at 
FUV wavelengths, in both cases for two different illumination angles /io. Note that the approximation used in deriving these curves becomes less accurate as yUo 
decreases, as the results start to develop an additional dependency on the viewing azimuth <p. (b) Effective albedos for diffuse transmission through a translucent 
foreground layer with optical depth 0.5 along the line of sight. Line types as in (a). 



approximation to these results is given by 

X(jt/, t) ^ 1 H- 0.75A, Y(^,T)^(l + 1.5A)e-^'^' (B14) 
where 

This approximation is accurate to < 10% for all the cases 
covered by Mayers (1962) (r - 0.1-5, mo = 0.5-1). To extend 
this result to forward-throwing phase functions, we assume 
that the anisotropy can be neglected for all orders of scattering 
higher than the first, so that equations (B14) and (B13) may be 
directly combined. 

Example results for scattering from foreground dust in this 
approximation are shown in Figure I2b, assuming (p - 45° and 
a Henyey-Greenstein form for the scattering phase function: 

1-9' 

(1 + r - ^g/^sV 

where 

iJs ^ ^l^io + {I - fi^y'Hi - fil)"^ cos 0. (B17) 

The results are shown for a fixed value of the optical depth 
measured along the line of sight, T/fi, since it is this quantity 
that is constrained by observations of the extinction in the 
neutral veil. Therefore, the actual thickness of the layer t goes 
to zero as yu — » and no limb brightening is seen. Instead, 
the albedo tends to have a maximum when fi ^ fio since this 
maximises O for the small values of (p considered here. It 
can be seen that the effective albedo is generally of order 
■wqt/ii, although it can be several times larger than this for 
favorable combinations of fi, //q, and (p that give sufficiently 
small scattering angles (/is > 0.8 for the optical-band grain 
properties, or /Vj > 0.9 for the FUV-band grain properties). 

B.3. Variations within the nebula of the illumination and 
viewing angles 

The effective scattering albedos derived in the previous sec- 
tions are strong functions of the angle of illumination of the 
scattering layer ji^ and of the observer's viewing angle //, with 
the albedo generally being highest when the illumination is 
close to face-on (/iq - 1) and the view is close to edge-on 



(ju ^ 0). It is therefore important to consider whether these 
angles vary systematically between the core and the outskirts 
of the nebula. This depends critically on the large-scale geom- 
etry of the scattering layers within the nebula. For instance, if 
the nebula were a simple hemispherical shell centered on the 
Trapezium stars (illustrated in Fig. 13a), then the illumination 
angle would be constant at ji^ - 1 while the viewing angle 
would vary from yu = 1 in the center to // = at the edge. On 
the other hand, if the nebula were a plane parallel layer (as in 
the models of Henney et al. 2005a and illustrated in Fig. \3b), 
then ji would be constant, whereas //q would vary from ^ 1 at 
the center to ^ towards the edges. In reality, neither of these 
simple geometries works well as model for the nebula. In par- 
ticular, the hemispherical-shell model would predict a constant 
ionized gas density and a surface brightness that increases with 
radius, both in violent disagreement with observations. On the 
other hand, the plane-layer model fails to explain the fine-scale 
structure see in many emission lines (e.g., O'Dell & Yusef- 
Zadeh 2000; Garcia-Diaz & Henney 2007), as well as the sharp 
edge of the EON. For an observational aperture that is larger 
than the angular size of the individual emission structures, the 
observed emission will be biased towards face-on illumination 
angles and edge-on viewing angles, simply because those are 
the cases that give the highest effective albedo. 



B.4. Breakdown of the infinite plane-parallel layer 
approximation 

The results of the previous sections assume that the scatter- 
ing occurs in a single plane-parallel layer of infinite lateral 
extent, which is a good approximation so long as the thickness 
of the each scattering layer and the displacements between 
them are much smaller than either the distance from the il- 
luminating source, or radius of curvature of the layer, or the 
size of the observational aperture. Obviously, these conditions 
will be violated to a greater or lesser extent in a real nebula, 
which will lead to a variety of aditional effects on the line 
ratios. The most important of these can be characterized as 
(1) ionization stratification, (2) differential pre- attenuation, or 
(3) limb-brightening limiting. We now discuss these in turn 
and show that none of them is likely to have an important effect 
on the observational results discussed in this paper. 
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Figure 13. Three simple models for the geometry of the nebula, showing 
how the illumination angle /^o and viewing angle /j vary with position in 
the nebula. In each case, the observer is located off the page to the bottom. 
(a) A hemispherical shell, {b) A nearly plane layer (c) An in'egular nebula 
consisting of many globule-like and bar-like features. 

Ionization stratification is the angular separation on the plane 
of the sky of the different scattering layers, such as the sep- 
aration of the Hy6 emission, which arises in the ionized gas, 
from the dust-scattered optical continuum, which arises pre- 
dominantly in the neutral PDR. This stratification is not visible 
in a true plane-parallel geometry unless the viewing angle is 
strictly edge-on (yu - 0), but for a finite geometry it will occur 
for < z/R where z is the separation between the layers and 
R is the smaller of the radius of curvature or the lateral extent 
of the layers. Although ionization stratification will produce 
fine-scale variations in the line ratios and equivalent widths, it 



will not affect the values given in Figure 1 unless the angular 
size corresponding to the inter-layer separation z is larger than 
the size of the observational sample regions. The sizes of the 
sample regions are listed in Appendix A of O'Dell & Harris 
(2010) and range from about 1 to 7 arcminutes, which are 
comfortably larger than the observed inter-layer separations 
in the regions within 7' of the Trapezium that are included in 
Figure 1 of this paper. 

Pre-attenuation is the reduction of the flux Fq incident on 
the scattering layer due to absorptions in material at smaller 
radii that does not contribute to the observed scattered intensity 
/. In the brightest regions of the nebula it can be shown that 
the majority of the emission in ionized lines such as H/S arises 
in a relatively thin layer near the ionization front (e.g.. Wen & 
O'Dell 1995), but there is also a more extended diffuse compo- 
nent to the emission, which becomes relatively more important 
at greater distances (Baldwin et al. 1991; Henney et al. 2005a). 
This pre-attenuation will afi'ect the line ratios and equivalent 
widths only if it is differential, that is, affecting one scattering 
process more than another For incident radiation in the optical 
and far-ultraviolet bands, the dominant absorption process is 
always due to dust, whereas for ionizing extreme-ultraviolet 
radiation it may be dust or hydrogen, depending on the local 
ionization parameter. In § B.4.1 below it is shown that, in the 
difi'use ionized gas responsible for the pre-attenuation, dust is 
the dominant opacity source in the EUV band also. Since the 
dust absorption cross section is very similar at optical and EUV 
wavelengths (Fig. 19 of Baldwin et al. 1991), pre-attenuation 
will have almost no effect on EW(Hy6), which is sensitive to the 
EUV/optical flux ratio. The dust absorption cross section in 
the EUV band is about 20-50% higher than in the visual band, 
so that pre-attenuation may affect EW([Ni],Corr), which is 
sensitive to the FUV/optical flux ratio. However, the total 
continuum optical depth to the [Ni]-scattering layer is only of 
order unity (see Fig. 3) and the optical depth of any diffuse 
pre-attenuating gas must be substantially less than this, so the 
effect is likely to be small. 

Limb-brightening is the increase in intensity of the emergent 
intensity as the viewing angle becomes more closely edge-on, 
due to the increased optical path through the scattering layer. 
In a strict plane-parallel approximation, the limb brightening 
does not saturate until the scattering layer is optically thick 
to the emergent radiation along the viewing direction, but any 
curvature in the layer will impose an additional limit on the 
degree of limb-brightening. This arises since the maximum 
path length through the layer is approximately 2 y/lRh, where 
R is the radius of curvature and h is the layer thickness, mean- 
ing that the maximum boost that limb-brightening can give the 
emergent intensity with respect to the face-on (jj - 0) value 
is of order 2 y/lR/h. Typical values of R/h vary from ^3.5 
for the HyS-scattering layer to ^ 100 for the [Ni] -scattering 
layer, giving maximum boost factors of ^ 5 and ^ 30, respec- 
tively. Since optical depth effects also limit the boost factor 
to a maximum of about 5 (see §§ B.2.1 to B.2.2 above), the 
extra limiting of limb-brightening by curvature effects will be 
unimportant, except arguably for H/3. 

B.4. 1. Relative importance of dust versus iiydrogen opacity at EUV 

wavelengths 

By using the equation of local photoionization equilibrium 
to rewrite the hydrogen photoabsorption rate in terms of the 
recombination rate, it is straightforward to show that dust will 
dominate the EUV opacity in ionized gas for densities less than 
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n' — FQO-,,„^JaB, where crj„„ is the EUV dust absorption cross- 
section. Taking crj„,, = 5 x lO^^^cm"^ H"^ (Fig. 19 of Baldwin 
et al. 1991) and using the ionizing luminosity of the Trapezium 
stars listed in Table 1, one finds n' ^ cm"-', where 

D is the projected distance from the Trapezium in arcminutes 
(assumed to be on average V3/2 times smaller than the true 
distance). Coincidentally, this equation for «' is very close to 
the reference line drawn on Fig. 6 of O'Dell & Harris (2010), 
which shows observationally derived electron densities as a 
function of distance, and from which it can be seen that n < 
n' for D < 2' but that n ~ n' for D = 2'-7'. The diffuse 
ionized gas in the interior of the H ii region is likely to have 
somewhat lower density than the mean densities derived from 
line ratios. Therefore, we conclude that dust is the dominant 
opacity source for EUV radiation in the difiiise ionized gas at 
all radii covered by our observations. 

Note, however, that this does not mean that dust is the domi- 
nant EUV opacity source in the H ii region as a whole. In fact, 
only 10-20% of the ionizing photons are absorbed by dust, but 
the hydrogen absorption is weighted towards the edge of the 
H II region, rather than the diffuse interior gas. 

C. NON-THERMAL LINE BROADENING OF [N i] AND 
OTHER LINES 

In the Orion Nebula, as elsewhere in the interstellar medium, 
significant non-thermal linewidths are observed to be ubiqui- 
tous in all gaseous phases (O'Dell 2001; O'Dell et al. 2003). 
This is shown in Table 12, which collates measurements from 
the literature of linewidths 6V and gas temperature T for vari- 
ous emission lines in the central Orion Nebula, ranging from 
fully molecular to fully ionized species. The expected ther- 
mal FWHM is SVth = 0.214 ^/(TJA) km s"' where T is the 
temperature in K and A is the atomic weight of the emitting 
species in units of the proton mass nip. This is subtracted in 
quadrature from the total width to give the non-thermal broad- 
ening component: <5y„th = {SV^ - ^V^^^y^^. It can be seen that 
the non-thermal component dominates over the thermal in all 
cases and increases in magnitude from about 3 km s"' in fully 
molecular gas up to about 15 km s ' in the fully ionized gas. 
If the non-thermal broadening is truly due to gas motions, then 
an approximate characteristic Mach number of these motions 
can be calculated as M = 0.5 5Vnth/Cs, where - ^{kT Ijinip) 
is the sound speed and ji is the mean mass per particle. This 
Mach number is shown in the last column of the Table, and in 
contrast to the linewidth it decreases with increasing ioniza- 
tion of the gas: the non-thermal motions are highly supersonic 
in fully molecular gas, slightly supersonic in the PDR, and 
slightly subsonic in the ionized gas. 

Optical and infrared emission lines from H ii regions and 
PDRs, such as the majority of those listed in Table 12 are usu- 
ally optically thin. Therefore, the observed line widths give no 
information about the spatial scales at which the broadening 
mechanism operates. On the other hand, for optically thick 
lines, such as the EUV lines that are responsible for pumping 
the [Ni] emission, one can divide potential broadening mech- 
anisms into two categories: microscopic and macroscopic, 
according to whether they occur at scales that are smaller than 
or larger than the relevant photon mean free path. Of the two, 
only microscopic mechanisms act to broaden the absorption 
profile and so affect the radiative transfer of the line, whereas 
macroscopic mechanisms simply act to broaden the emergent 
intensity profile. 



The line broadening that we derive for the pumping lines 
in order to explain the observed optical [N i] line brightness 
(see Figure 5) is similar to, but smaller than, the broadening 
observed in the lines themselves (Table 12), implying that the 
non-thermal broadening mechanism must be microscopic in 
nature. In other words, it should occur on scales of less than 
10'** cm, which is the approximate mean free path of the 954 A 
pumping line. 

In the H ii region, transonic turbulence is expected to be 
driven at the scales of photoevaporation flows from dense glob- 
ules and filaments (Mellema et al. 2006; Arthur et al. 2011; 
Ercolano et al. 2011). The most vigorous photoevaporation 
flows only occur at scales larger than about 10% of the Hii 
region radius (Henney 2003). In the Orion nebula, the closest 
approach of the ionization front to the ionizing stars is about 
0.2 parsec (Wen & O'Dell 1995), so we assume that turbu- 
lent velocities of amplitude 1 1 km s"' are present at scales of 
0.02 parsec, or 6 x lO'^ cm. In a Kolmogorov-type turbulent 
energy cascade, velocity differences scale with separation I 
as 6v ~ {^^^. Therefore, on the scale of the Ni mean free 
path the turbulent broadening should be only approximately 
1 km s"', which is much smaller than our derived microscopic 
broadening, which means that this is not the mechanism we 
seek. 

A further potential source of broadening is the systematic 
acceleration of the gas as it is dissociated, heated, and ionized 
by the advancing front. However, the [N i] emission arises in 
regions where the gas is still predominantly neutral, with tem- 
perature ranging from 1000 to 5000 K (see Fig. 3), for which 
the velocity increase is expected to be small since the greater 
part of the gas acceleration occurs in the warmer, partially 
ionized zone where the [Oi] lines arise. For example, a plane- 
parallel model of a D-critical ionization front (Eqs. (A5)-(A8) 
of Henney et al. 2005b) implies a total broadening FWHM for 
the [Ni] lines of less than 2 km s"' by this process. 

We therefore see that none of the broadening mechanisms 
that have been succesfully invoked to explain the observed 
widths of neutral and ionized coUisional lines are successful 
in explaining the observed characteristics of the fluorescent 
[N i] lines. Other potential mechanisms such as instabilities of 
the ionization front itself (Williams 2002; Whalen & Norman 
2008) are not promising either, since they are unlikely to pro- 
duce microturbulence at a sufficiently small scale. Neither can 
broadening due to dust scattering in the neutral veil (Henney 
1998) be the explanation, since this would have no affect on 
the radiative transfer of the N i pumping lines. 

A more promising mechanism for generating turbulent veloc- 
ities on a very small scale is the action of thermal instabilities 
(Koyama & Inutsuka 2002) in the shocked neutral layer that 
precedes the ionization front. It is suggestive that the tem- 
perature range over which the [Ni] lines form in the PDR 
(1000-3000 K) is similar to the range over which the ISM is 
known to be thermally unstable (Field et al. 1969). If such 
an instability were to occur in the PDR, then the smallest 
scale at which fragments could arise (and hence turbulence 
be driven) is given by the Field length Ap (Field 1965), which 
represents the scale below which temperature fluctuations will 
be smoothed out by thermal conduction. Assuming saturated 
conduction by free electrons (Zel'Dovich & Raizer 1967), one 
finds a value of /Ip - (10'^/«) cm, which is roughly 100 times 
smaller than the thickness of the [N i] pumping layer The role 
of thermal instability in generating the required microscopic 
non-thermal broadening therefore merits further investigation. 
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Table 12 

Non-thermal line widths of different gas phases in Orion 

Atomic Mean Gas Sound Line width FWHM Mach 

Weight mass temperature speed Total Thermal Non-thermal number 

Species A j^i T (K) <;v(kms"'j ()V'(kms"') ()V'ih(kms"') (>"V'„,i, (km s"' ) M 

~CO 28 236 40 ± 10 0.37 ± 0.05 3.0 ± 0.5 0.26 ± 0.03 2.99 ± 0.50 4.0 ± 0.9 

[Cii] 12 1.30 400 ± 100 1.59 ± 0.20 4.0 ±1.0 1.23 ±0.15 3.81 ± 1.05 1.2 ±0.4 

[Ni] 14 1.30 2000 ±500 3.56 ± 0.45 6.0 ±3.0 2.56 ± 0.32 5.43 ± 3.32 0.8 ±0.5 

[Oi] 16 0.90 9500 ±500 9.33 ± 0.25 12.6 ±2.4 5.21 ±0.14 11.47 ±2.64 0.6 ±0.1 

[Nn] 14 0.68 9000 ± 500 10.45 ± 0.29 16.3 ± 3.5 5.42 ±0.15 15.37 + 3.71 0.7 + 0.2 

[Qui] 16 0.65 8400 ± 500 10.33 ±0.31 15.5 ± 4.8 4.90 ±0.15 14.71 ± 5.06 0.7 ± 0.2 
References. — CO: Wilson et al. (2001); [Cii]: Boreiko et al. (1988); other Unes: Baldwin et al. (2000); Garci'a-Di'az et al. (2008) 
Note. — All linewidths have been corrected for instrumental broadening. 



